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Abstract 

In quantum physics, recent investigations deal with the so-called "quantum trajec- 
tory" theory. Heuristic rules are usually used to give rise to "stochastic Schrodinger 
equations" which are stochastic differential equations of non-usual type describing the 
physical models. These equations pose tedious problems in terms of mathematical 
and physical justifications: notion of solution, existence, uniqueness... 

In this article, we concentrate on a particular case: the Poisson case. Random 
Measure theory is used in order to give rigorous sense to such equations. We prove 
existence and uniqueness of a solution for the associated stochastic equation. Fur- 
thermore, the stochastic model is physically justified by proving that the solution 
can be obtained as a limit of a concrete discrete time physical model. 



Introduction 



Many recent developments in quantum mechanics deal with "Stochastic Schrodinger Equa- 
tions". These equations are classical stochastic differential equations (also called "Belavkin 
equations") which describe random phenomena in continuous measurement theory. The 
solutions of these equations are called "quantum trajectories", they give account of the 
time evolution of an open quantum system undergoing continuous measurement. 

Usually, in Quantum Optics or Quantum Communication, indirect measurement is 
performed in order to avoid phenomena like Zenon effect. The physical setup is the one 
of a small system (an open system) interacting with an environment, the measurement is 
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then performed on the environment. In the literature, stochastic Schrodinger equations are 
expressed as perturbations of the Master Equation which describes normally the evolution 
of the small system without measurement. 

Belavkin equations are stochastic Schrodinger equations describing the evolution of 
qubit, i.e two level system. Essentially there exists two type of equations. 

• A diffusive equation 

d Pt = L(p t )dt + [p t C* + Cp t - Tr [p t (C + C*)] pt] dW t (1) 
where W t designs a one-dimensional Brownian motion 

• A jump-equation 

J(Pt) 



dp t = L(p t )dt 



Tr[j(pt 



Pt 



(dN t -Tr[j(p t )]dt) (2) 



where N t is a counting process with intensity J" * Tr [(J'(ps)] ds (the different operators 
are described in the paper). 

In this article, we shall focus on the case of jump-equation ( |29|) . the diffusive case ([I]) is 
treated in details in [22] . The equation (2) pose tedious problems in terms of mathematical 
justification. In the literature, the question of existence and uniqueness of a solution 
is not really treated in details. Actually classical theorems cannot be applied directly. 
Furthermore, the way of writing and presenting this equation is not clear. Indeed, how 
can we consider a driving process which depends on the solution? With the expression (2), 
there is no intrinsic existence for such process. Even the notion of solution is then not clear 
and a clearly probability space must be defined to give rigorous sense for such equations. 

Regarding the physical justification of Belavkin equation model, heuristic rules are 
usually used to derive these equations (see [§]). In order to obtain them rigorously, often 
a heavy analytic machinery is used (Von-Neumann algebra, conditional expectation in 
operator algebra, Fock space, quantum filtering...). This high technology contrasts with 
the intuition of heuristics rules. In this article for the very first time, the continuous model 
(|29|) is justified as a limit of a concrete discrete model of quantum repeated measurements. 
This approach, which has been already developed for the diffusive model ([T]) in [22], is 
based on the model of quantum repeated interactions (see [3] for all details). The setup 
is the one of a small system interacting with an infinite chain of quantum system. Each 
piece of the chain interacts with the small system, one after the others, during a time h. 
After each interactions a measurement is performed on the piece of the chain which has 
interacted. The sequence of measurements induces then random perturbations of the small 
system. This is described by a Markov chain depending on the parameter h. Such Markov 
chain is called a discrete quantum trajectory. In this article, by using a renormalization for 
the interaction (based on the result of |3J) and by considering the limit h goes to zero, we 
show that the quantum trajectory corresponding of the jump-equation can be obtained as 
continuous limit of particular discrete quantum trajectories. 
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As it is mentioned, similar results concerning existence, uniqueness and approximation 
for the diffusive equation is expressed in [22J. The particularity of the jump-equation 
concerns the counting process (N t ). In the diffusive case, the existence and uniqueness 
of solution concerns a strong notion of solution. In the jump context, in a first time, 
we use the notion of Random Poisson Measure to give a sense to the equation and in 
a second time, we deal with existence and uniqueness of a weak solutions (i.e path by 
path). Concerning the convergence result in the diffusive case, the arguments are based on 
a convergence result for stochastic integrals due to Kurtz and Protter ([IE], In the 

jump case, such techniques cannot be applied (essentially due to the stochastic intensity, 
see section 3). The result is then based on Random Coupling Limit method which allows 
to compare discrete quantum trajectories and continuous quantum trajectories. Finally 
to conclude to the convergence, we need the intermediary result of convergence of Euler 
scheme for the jump-equation. This result is also shown in this article because classical 
results cannot be applied in this context. 

This article is structured as follow, we can Section 1 is devoted to remind the discrete 
model of repeated quantum measurements based on quantum repeated interactions models 
and indirect measurement. Next, we make precise the model of a qubit in contact with a 
spin chain; this model corresponds to the effective discrete model of Belavkin equations. 

In Section 2, we investigate the study of continuous model of jump-equation. We present 
a way to define a concrete probabilistic framework using random Poisson measure theory. 
Next, we solve problems of existence and uniqueness. The way, we describe the solution, 
is next used in Section 3 to prove the convergence of discrete model. 

In Section 3, we prove the convergence of a special discrete model to the jump continuous 
model. The arguments, we use in the jump case, are totally different of the diffusive case 
[22] . Indeed the one, used in [22], are based on results of Kurtz and Protter ([18]. [T9j ) and 
cannot be applied in this context. In order to prove the convergence result, we use explicit 
comparison by using random coupling method. Furthermore, we need to show another 
result of approximation of jump-equation: the Euler scheme convergence. This result is 
next use to prove the final convergence result. 

1 Quantum repeated measurements: A Markov chain 
1.1 Quantum repeated measurements 

In this section, we present the mathematical model describing the quantum repeated mea- 
surements setup. Such model is treated in details in [22]- We do not remind all the details, 
we just present the rules and the Markov property of discrete quantum trajectories. 

The model is based on quantum repeated interactions model. As it is mentioned in the 
introduction, we consider a small system Tio in contact with an infinite chain of quantum 
systems describing the environment. All the pieces of the chain are identical and indepen- 
dent of each other; they are denoted by TC. One after the others, each copy of 7i interacts 
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with Hq during a time h. After each interaction, a measurement is performed on the copy 
H, it involves then a random perturbation of the state of Ho- 

Let start by describing a single interaction between Ho and H and the indirect mea- 
surement on H. Let Hq and H be finite dimensional Hilbert space. Each Hilbert space is 
endowed with a positive operator of trace one. Such operators represent states of quantum 
systems. Let p be the initial state on H and let (3 be the state of H. 

The coupling system is described by H <E> H. The interaction is described by a total 
Hamiltonian H tot which is a self adjoint operator defined by 

H tot = H ®I + I®H + H mt . 

The operators H and H are the Hamiltonians of Ho and H; they represent the free 
evolution of each system. The operator ifj nt is called the interaction Hamiltonian and 
describes the energy exchanges between the two systems. The Hamiltonian H tot gives rise 
to a unitary operator of evolution 

U = e ihHtot , 

where h is the time of interaction. Hence, after the interaction, in the Schrddinger picture, 
the initial state p <g> (3 on Ho <g> H becomes 

p = U{p®(3) U*. 

Let us now describe the measurement of an observable A of H. An observable is a self- 
adjoint operator and we consider its spectral decomposition 

i=0 

Actually, we consider the observable / <8> A on H$ <8> H. According to the law of quantum 
mechanics, the measurement of / <g) A gives a random result concerning the eigenvalues of 
I <g) A. It obeys to the following probability law 

P[to observe A,] = Tr[pl ® P { ]. 

After the measurement, if we have observed the eigenvalue Aj the reference state p becomes 

... I®PipI®Pi 

Pl{l) = r 1 ■ 

Tr [ p I ® Pi \ 

Such phenomena is called Wave Packet Reduction Principle. The state pi(i) is then the 
new reference state of Ho <8> H conditionally to the observation of Aj. 

In general, one is only interested in the evolution of the small system Ho. We use the 
partial trace operation to define a state on Ho from a state on Ho ® H. The partial trace 
operation is defined as follows. 
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Definition-Theorem 1 Let a be a state on the tensor product Ho®H. Hence there exists 
a unique state on Ho, denoted E [a] ; which is characterized by the property 

VX e B(H ) Tr Ho i V X] = Tr no ^ H [a(X ® I) ], 

where Tr-^ corresponds to the trace of operator on Hq and similar definition for Tr-^®?^. 

Now, let define the state p\(.) on Ho by 

Px{i) = E [//i(i)] , i = 0,...,p. 

The state Px(-) is then the new reference state of Ho describing the result of one interac- 
tion and one measurement. It represents a random state; each state pi(i) appears with 
probability pi = P[to observe A*]. The probability pi represents then the probability of 
transition from the state p to the state pi(i). 

As a consequence, Ho is now endowed with the state p\ and a second copy of H can 
interact with Ho- In the same way, a measurement of I <8> A is then performed and by 
taking the partial trace, we get a new random variable p 2 with similar transition probabil- 
ities (from pi to p 2 ) an d so on. Hence we can define a random sequence (pf.) of states on 
Ho ; this sequence is called a discrete quantum trajectory. It describes the random evolu- 
tion of the state of Ho undergoing quantum repeated interactions and quantum repeated 
measurements. 

From the description of the rules of one interaction and one measurement and by 
construction of the random sequence (pk), the following proposition is straightforward. 

Proposition 1 The discrete quantum trajectory (p^ is a Markov chain. The transition 
are described as follows. If pk = Xk then p^+i takes one of the values 

_ \l®PiU(xk®P)U*I®P^\ . . 

En f ; ; t~ , 1 = U, . . . , V, 

lTr[U(xk®P)U*I®Pi] J' 
with probability Tr [ U(xk ® /3)U* Pi ] . 

A complete justification of this proposition and a construction of an appropriate prob- 
ability space for the Markov chain can be found in [22] • All the details are not necessary 
in the following. 

The next section is devoted to the study of a special case of a qubit in contact with a 
chain of spin. 

1.2 A qubit interacting with a chain of spin 

The mathematical setup describing such model is represented by Ho = H = C 2 . In this 
section, we make explicit a discrete version of Belavkin equations. 

Let us start by describing an evolution equation for the state (pk) in this context. In the 
case of C 2 , interesting observable owns two different eigenvalues, that is A = XqPq + A]Pi. 
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The description of transitions of the Markov chain can be expressed by the following 
equation 

_ £o(pfc) lfc+ i , £i(pfc) -.fc+i /„x 

Pk+l ~ Tr[£ (p*)] ° + Tr[C 1 (p k V l1 ^ ( " >j 



where for % G {0, 1} terms Ci(p k ) corresponds to the "non-normalized" transition of Pk+x, 
that is Ci(pk) = E o [l®P i U(p k ®0)U*I®P i ]. We denote Pk+1 = Tr[C (p k )] =l-q k+1 . 
Moreover, the term 1q +1 corresponds to the random variable which takes the value 1 with 
probability p k+ i and with probability q k +i] it corresponds to the observation of Ao at the 
k+l-th measurement. The inverse holds for the random variable 

In order to make more precise the evolution of the state (p k ), we have to express in 
a more explicit way the terms Ci(p k ). For this, we introduce an appropriate basis. Let 
(Q,X) denote an orthonormal basis of C 2 . For 7^ ® Ti, we consider the following basis 
Q (gi Q, X (gi Q, Q (gi X, X (gi X. In this basis, the unitary-operator U can be written as 



U 



Loo Lqi 
Lio L\\ 



where each are operators on 7i . For the state j3, let choose the one dimensional 
projector on Q, that is 

Furthermore if for i G {0,1}, the projector Pj = {p kl ) ky i=o,i in the basis ($7, X), after 
computation we get 

d(p k ) = Eo[l®PiU(p k ®P)U* I®Pi] 

= Poo Loo Pk Lq + poi -^oo Pk L* + p\ Lio Pk Lq + p l n L 10 p k L\ (4) 

In order to compare the discrete evolution and the continuous evolution in Section 3, 
we modify the way of writing the equation ([3]) by introducing new random variables. For 
all k > 0, let put u k+ i = and define the random variable 

v _ v k+ i — q k+ i 
Afc+i — — 

y/Qk+iPk+i 

We define the associated filtration on E N 

T h = a(Xi,i < k). 

So by construction we have E[X k+ i/ J 7 ^ = and E[X| +1 /jF fc ] = 1. The random variables 
(X k ) are then normalized and centered. In terms of {X k ), the discrete evolution equation 
for the discrete quantum trajectory becomes 



p k+ i = C (p k ) + Ci(pk) + 



Qk+i r , \ | Pk+l r i > 



Pk+l V Q k + l 



X k+1 . (5) 
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This way, the evolution of the discrete quantum trajectory appears as a random perturba- 
tion of the equation p k+ i = C (pk) + £i(pk) which describes actually the evolution without 
measurement. In [3], it is shown that + C\ is a completely positive application and that 
it is a discrete Master equation. Equations of type (J5j), obtained for different observable, 
can be then compared with Belavkin equations in terms of perturbation of Lindblad Master 
equations (see Section 2). 

Remark: It is worth noticing that such equations depends only on the expression of 
eigen-projectors and do not depends on the value of eigenvalues. 

The following section is devoted to the study of the continuous Belavkin equations with 
jump. 



2 The jump Belavkin equation 

First rigorous results and mathematical models describing evolution of systems undergoing 
quantum measurements are due to Davies in [8]. From his works, one can derive heuris- 
tically the Belavkin equations. A heavy background is necessary in order to obtain these 
equations in a rigorous way (quantum filtering theory |3], instrumental process |27j). In 
this article, the continuous stochastic model is rigorously justified as a limit of the discrete 
process given by the equation ([5]). 



In this article we focus on the jump-equation 



dp t = L(p t )dt + 



J(Pt) 



Tr[j(p t 



Pt 



dN t -Tr[j(p t )]dt 



(6) 



where N t is assumed to be a counting process with intensity J Q Tr\j(p s )\ds. 

Before to study this equation, let us speak briefly about the different operators appear- 
ing in the expression. In quantum physics, the operator L is called the Lindbladian of the 
system. This is a classical generator of the dynamic of open quantum systems. It gives 
rise to the Master equation (see the remark at the end of Section 1) 



j t p t = L( Pt ) = -i[H, - i {C*C, p t } + C Pt C* 



(7) 



where C is any 2x2 matrix, the operator Hq is the Hamiltonian of the qubit. In the 
equation ([6]), the operator J is defined as J(p) = C pC*. 



2.1 The probability framework of the jump- equation 

The way of writing equation ([6]) and defining the process (N t ) is not absolutely rigorous in 
mathematical point of view. Indeed in the definition, the process (N t ) depends on existence 
of the solution of the jump-equation. However to treat the problem of the existence of a 
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solution, we have to consider firstly the definition of the driving process of the stochastic 
equation. 

Actually, we cannot consider the existence of the process (N t ) without the existence of 
the solution (p t ) and reciprocally. It imposes the following definition of solution for the 
jump-equation. 

Definition 1 Let (Q, T, J-'t, P) be a filtered probabilistic space. A process-solution of © 
is a cadlag process (p t ) such that there exists a counting process (N t ) with predictable 
compensator (or stochastic intensity) 



Tr[J(p s J)\ds 

and such that the couple (pt, N t ) satisfies almost surely 



Pt = Po + 



o 



L{p B -)- J(p s -)+Tr[j(p s 



Ps 







ds + 


/'[ 

Jo . 



J(Ps 



Tr[J{p s J)\ 



Ps 



dN K 



This notion of solution comes from considerations of Jacod and Protter in [Hj [11] . 

After this definition, the next step consists in constructing a process iV with stochastic 
intensity. For this, we use the general theory of Random Measure (for all details, see [UJ 
or Let us introduce this notion. 



Definition 2 Given a filtered probability space (Q, T , J-'t, P), a random measure is a family 
of measure \i = (p(u, G fi) on (R + x R d , B(R+) ® B(K d )). 
A random measure is said to be integer valued if 

1. For all oo G p(uj,tx R d ) < 1. 

2. For all A e B(R+) ® B(R d ), the quantity p(A) is valued in N[J{+oo}. 

Definition 3 A random Poisson measure on (Q, J 7 , J-'t, P) is a integer valued measure that 
verifies 

1. The measure m(A) = E(p(A)) on B(R + ) <S> B(R d ) is non atomic. 

2. m(0 x R d ) = 0. 

3. If t G R+ and if A4 G B(]t, +oo[), i = 1, . . . ,1 are two by two disjoint with m(Ai) < 
+00, the random variables p(Ai) are mutually independent and independent from Tt- 

The measure m is called the intensity of the random Poisson measure p. 

The following theorem shows how the random measure theory is used to construct the 
process (Nt)- 
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Theorem 1 Let (Q, J 7 , Tt, P) be a filtered probability space which supports a random Pois- 
son measure p on R x R with intensity dt®dx. Every process-solution of the equation 



Pt 



Po + 



+ 



L(p s _)+Tr[j(p s _)] Ps --J(p 

J(Ps-) 



ds 



JR 



Tr[J(p s 



lo<x<Tr[J( Ps „)}H(ds, dx) 



(8) 



is a process-solution of equation @ satisfying Definition^ For the process (Nt), we put 



N t 



JR 



lo<x<Tr[J(p s -)]Kds, dx). 



(9) 



A general form of this theorem can be found in [13]. In this theorem, there are two 
parts. On the one hand, we must prove that the process given by © is well defined, that 
is, it is a non-explosive process. On the other hand, we must prove that any solution of 
equation (jSJ) satisfies Definition [TJ 

The non-explosive property of (Nt) is related to the boundness character of the stochas- 
tic intensity t — > Tr[J~(p t -)]. Here, a straightforward computation shows that there exists 
a constant K such that for all state p, we have < Tr[J(p)] < K. It implies directly that 
for all t, the quantity Tr[J(p t _)} = \\m s<ts ^ t Tr[J(p s )] satisfies 

< Tr[J(p t ^)] < K, 

the intensity is then bounded. This property appears also in the proof of the theorem. 
Proof: Let show that the counting process (Nt) given by ([9]) is non explosive. For all 
cadlag matricial process (X t ), we define the time explosion 



T 



x 



inf it : N 



x 



+oo 



Let show that, if (p t ) takes values in the set of states, the explosion time T p = oo almost 
surely. For this, we introduce the following stopping times 



T n = inf U ^ Tr[J(p s _)]ds > n\ 



n > 1. 



It was clear that J Q Tn Tr[J(p s J)\ds < n and thanks to the property of Poisson random 
measure 



E 



Nt 

J- n 



E 
E 



R 



l 0<x<Tr[J{p,J)\ 



dsdx 



Tr [J(p s -)}ds 
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As a consequence we have T n < T p almost surely. It is then sufficient to show that 
limT„ = oo. The property < Tr[J{p t )] < K implies §*Tr[J{p s Jj\ds < oo, hence we 
have limn^ooTn = +00. Finally we have construct a counting process without explosion 
for all cadlag process which takes values in the set of states. Concerning the property that 
a process-solution of (jHJ) satisfies Definition (TJ this result follows from the construction of 
(N t ). □ 

Remark: This remark concerns a way to construct a random Poisson measure with uni- 
form intensity measure. Let (Q, P) be the probability space of a Poisson point process 
iV on R x R. The natural random Poisson measure attached with iV is defined for all 
A e B(R) <8> B(R) by 

p(.,A) = N(.,A). 

For all Borel subset A G B(R) <g> B(R), we have E[N(.,A)] = X(A) where A denotes the 
Lebesgue measure. The Poisson random measure p satisfies then the conditions of Theorem 
[TJ This particular random Poisson measure is used in Section 3 to realize continuous 
quantum trajectories and discrete quantum trajectories in a same probability space and to 
prove the convergence result. 

After defining this clear probability framework, it remains to deal with problems of 
existence and uniqueness of solution. 



2.2 Existence and uniqueness 

Let (Q, J 7 , Tt-, P) which supports a random Poisson measure p. In order to solve problems 
of existence and uniqueness of solution for equations of type (jHJ), a way consists firstly in 
solving the ordinary differential part 

Pt = p + [ \L(p s _) + Tr[j( Ps _)]p s _-J(p s _)}ds, (10) 
Jo 1 J 

and secondly defining times of jump. Classical theorems are usually based on Lipschitz 
property (cf [20], [H] or [7]). However in our case, the ordinary differential equation 
part contains non-Lipschitz functions. Hence the application of classical theorems is not 
st raight forward . 

The first step is then solving a non-Lipschitz ordinary differential equation (1101) . For 
this, we use the fact that this equation preserves the property of being a pure state. Remind 
that pure states are a particular class of states which are one dimensional projectors. The 
property of preserving pure states for (ITU]) is expressed as follows. 



Proposition 2 Let x be any vector of norm one in C 2 , if the Cauchy problem 

dx t = [—iH - \C*C + \r]t\ x t dt , , 

Xq = x 

where rj t = (x t , C*Cx t ) has a solution then \\x t \\ = 1 for all t > 0. 
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Furthermore the process (p t ) of one- dimensional projector defined by p t = P{ Xt } f° r a ^ 
t > is solution of Cauchy problem 



dpt 



L(p t )+Tr[j(p t )]p t -J(pt 



dt 



(12) 



Proof: Let (x t ) be the solution of fill I) , thanks to the fact that H is self-adjoint and 
r\t = (x t ,C*Cx t ) = (Cxt,Cx t ), a straightforward computation gives 



d_ 

If 



(x t ,x t ) 



(j t x t ,x t ) + (x t ,±x t ) 



( 



\X t , 



-in, - hc*c + \ m 



x t , x t ) + 

Xt) 



1 1 / X 

-lj?k + -^lt\xt,x t ). 



As a consequence, if (x ,x ) = 1, then (x t ,x t ) = 1 for all t. Let p t = P{ Xt }, f° r an V, we 
have 

pt y = (x t ,y)x t . 

We can derive dpty and we have 



d 
dl 



PtV 



{^xt,y)x t + (x t ,y) — x t 



( 



(x t , 



Pt 



-iH Q - l -C*C + ^ 

-iH - ^crc + \m 



x t ,y)x t + (x f ,2/) 
y)x t + 



-iH - \c*C + ~ m 



-iH - l -C*C + 



x t 

{x t ,y)x t 



y + 



-iHo ~ \c*C + \vt 



PtV 



= L(p t )+Tr[j(p t )]p t -J(p t 
and the result follows. 



□ 



With this expression in term of pure states, we can show that equation ( flOl) admits a 
solution. For this, we need the following characterization of pure states in C 2 . 

Lemma 1 Let p be a state on C 2 . // there exists a vector x G C 2 such that (x; px) = 0, 
the state p is a one dimensional projector. 

We do not give the proof, it is a simple linear algebra result. 

The following proposition expresses the result concerning the existence of a solution for 

(HDD. 



n 



Proposition 3 Let p be any state, the Cauchy problem 



dpt 



L(p t )+Tr[j(p t )]p t -J{p t ) alt 



(13) 



P 



has a unique solution defined for all time t. 

Furthermore, if there exists to such that pt is a one dimensional projector, the solution 
of Q3l) after t is valued in the set of pure states. 

Proof: As the coefficients are not Lipschitz, the Theorem of Cauchy Lipschitz cannot be 
applied directly. However the coefficients are C°°, so locally Lipschitz and we can use a 
truncation method. The ordinary equation is of the following form 



where / is C°° and f(A) = L(A) + Tr[J(A))A - J (A). We define the truncation function 
ip from R to R defined by 



For a matrix A = (%•), we define by extension pk(A) = <p fc (i£e(ciy)) + i<fj c (Im(a i j)). Thus 
/ o (p k is Lipschitz. Now we consider the truncated equation 



The Cauchy-Lipschitz Theorem can be applied because f o(p k is Lipschitz and there exists 
a unique solution t i— »■ p k) t defined for all t. 
Now, we can define the following time 



As po is a state, if k is chosen large enough we have T k > and for all t < Tk, we 
have <fk{pk,t) — Pk,t- Thus t \-+ pk± is the unique solution of the ordinary equation (TTUT) 
(without truncation) on [0,7V]. The classical method in order to solve an equation with 
non Lipschitz coefficients is to put T = lim^ T k and to show that T = oo. 

Here the situation is more simply because if po is a state, we can show that the solution 
is valued on the set of states. Indeed, as ||p|| < 1 when p is a state, we have for example 
<f2(p) = P- Let show that the solution obtained by truncation is valued in the set of states, 
it will imply that T 2 = oo. We have to show that the solution is self adjoint, positive and 
of trace one. 

Let start with the self adjoint property and trace property. On [0, T 2 ], as the ordinary 
differential equation is Lipschitz, we can solve it by Picard method. Let define 



dpt = f(pt)dt 




-k if x < —k 
x if —A; < x < k 
k if —k < x < k. 



dpk,t = f ° <Pk(pk,t)dt. 



T k = inf{t,3(ij)/ \(Re(aij(pk,t)))\ = k or \(Im(aij(p k ,t)))\ = k}. 




P«(0) + J Q f o (fi k (pn(s))ds 
P- 



(14) 
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It is easy to show with the right definition of / that this sequence is self adjoint with trace 
one. Hence it proves (by convergence of the Picard method) that for all t < T 2 , the matrix 
P2,t is self adjoint with trace one. 

It remains to prove the positivity property which cannot actually be proved by the 
previous Picard method approach. This condition is nevertheless a consequence of Propo- 
sition [5] and Lemma [TJ We must prove that for all y £ R 2 and for all t < T 2 , we have 
(y,P2,ty) > 0. Let define 

T° = inf {t < T 2 , 3y £ R 2 /(y, p %t y) = 0}. 

If T° = T 2 , an argument of continuity for solution of ordinary differential equation 
implies that (y, pifij) > for all t < T 2 and all y £ C 2 . The solution is then valued on the 
set of states. 

If T° < T 2 , by continuity there exists x £ M 2 such that (x,p2 t T°%) = and for all 
t < T° and for all y £ R 2 , we have (y,P2,ty) > 0. It means that on [0, T°], the solution 
t | — ► P2,t is valued on the set of states. Moreover for some x, we have (x,p^ T0 ) = 0. 
Thanks to Lemma [TJ the operator p 2j T° is a one dimensional projector. Hence, we can 
now consider the ordinary differential equation with initial state pxo = po- We are face to 
the Cauchy problem (TTTi) which is equivalent to the problem ( fl3l) thanks to Proposition [21 
This problem can be also solved by truncation method, the fact that the norm is conserved 
implies that the solution is defined for all t (the truncation is actually not necessary). 
Thanks to Proposition [21 we have a solution for the initial Cauchy problem f[TH|) which 
defines a state-process. We have then finally proved that on [T°, T 2 ] the solution is valued 
on the set of states. 

On T 2 , the solution is then a state. A local argument and the uniqueness in the Cauchy- 
Lipschitz theorem allows us to conclude that T 2 = oo and that there exists a unique solution 
of the ordinary differential equation (1101) . □ 

This above proposition is essential in the proof of the following theorem concerning 
existence and uniqueness of the jump-equation (jHJ). 

Theorem 2 Let (fl, JF, jF t , P) be a probability space which supports a Poisson random 
measure p whose the intensity measure is dx <S> dt. Let po be any state, the jump-Belavkin 
equation 



Pt 

ft 



p + / U(p s _)+Tr[j r (p s _)]p s _- J(p B _) ds 
Jo 1 

J(Ps-) 



+ 



R+ [Tr[J(p s _: 



Ps- 



lo<x<Tr[j( Ps -)]P(ds } dx) (15) 



admits a unique solution defined for all time. The process-solution (p t ) takes values on the 
set of states. 
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Proof: Such equation is solved paths by paths. As the initial condition p is a state, 
Proposition [3] states the existence and the uniqueness for the Cauchy problem (|T3|) . Let 
define the first time of jump, we put 



P(l)t 
2i 



Po + Jo f(p( 1 )s)ds 
inf{t,A>f (1) >0}. 



(16) 



The first equality expresses the fact that the Cauchy problem admits a solution. The 
second equality defines a random variable Ti which is the stopping time of the first jump. 
We have 

N* 1 \u) = f i(u,G(j>,T l ,0)) = l, 
where G(p, t, s) = {(u, y) G R 2 /t < u < s, < y < Tr[J(p u )]}. The quantity p(u, G(p, t, s)) 
represents the number of point under the curves t — > Tr[J(p t )]. 

Now if T\ = oo, the solution of the jump-equation is given by the solution of the 
ordinary differential equation; there are no jumps. 

If T\ < oo, let construct the second jump time and explain how the solution of the 
jump-equation is defined. On [0,Ti[, we put p t = p(l)t and at the jump time T 1; we 
implement the value of the jump, that is, 



Pn = Pn- + 



Pn- 



Tr[J(p Tl -)Y 



[Tr[J( PTl - 

This matrix is well defined because the property T\ < oo implies Tr[J{pT 1 -)} > 0. More- 
over, the matrix is a state. The ordinary differential equation ffTUj) can be solved again 
with this initial condition. We can then define 

P(2)t 
P(2)t 
T 2 



p(l)t on [0,Tx[ 
Pt, + j£ f(p(2) a )dt 



mi{t> T 1; iVf (2) >N^>} 



rP(i)- 



(17) 



The random variable T 2 is then the second jump-time. If T 2 < oo, we have Tr[J'(p(2)j' 2 _)] > 
and we continue in the same way. The solution of the jump-equation is then defined as 
the solution of the ordinary differential equation between the jump times and at jump 
times we put the value of jump. 

More precisely, we define a sequence of random times T n and a sequence of processes 



P(n)t 

PT n -! 



p(n - l) t on [0,T„_i[ 
PT n _, + jj n l f{p{n - l) s )ds 



(18) 



inf{t >T n - X ,Np n) >N£ n - L >} 



yp(n-l) - 



All the processes are well defined because the matrices defined at each jump time are 
states and the Cauchy problem can be solved. The sequence of random stopping times 
(T n ) satisfies T n+ i > T n on the set {T n < oo}. In a classical way, we put 



lim T n 
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and we can define the process-solution (p t ) of the jump-Belavkin equation on [0, T[. For 
all t < T, we put 

Pt = p(n) t on [0,T n [. (19) 

This process is clearly a solution of the jump-Belavkin equation ( fl5l) and it is valued on 
the set of states. The uniqueness is implied by the uniqueness of the solution of Cauchy 
problem (cf Proposition [3]). Moreover any other solution is forced to have the same random 
jump-time, it implies the uniqueness. 

In order to finish the proof, we must show T = oo a.s. This random time is the 
explosion-time of N p , we cannot directly apply the result of the theorem 3 because the 
definition for all t of (p t ) is not already proved. However (p t ) is a state valued process, so 
we have Tr[J'{p t )} < K for all t then 



E 



< E 



NS 



< Kn. 



Furthermore N^ pAn 



P[T < n] 



p on {Tp < n}, it follows that pP[T p < n] < Kn, then we have 
for all n and the result is proved. □ 



The proof of theorem gives an explicit way to construct the solution and the random 
times of jumps. This construction can be used in a more general context. In the literature, 
the problem of existence is sometimes taken into account by a way using filtering method. 
A starting point is using a linear equation and next a change of measure. However this 
method do not give an explicit construction of the solution whereas we need it to prove 
the final convergence result. 

In order to prove the final result, we use random coupling method. For this, we need an 
explicit realization of the process (p t ). For this, we use Poisson point process theory. Let 
consider the probability space (fl, T, J-'t, P) of a Poisson point process JV on R x R. Let 
denote N(uj, ds, dx), the differential increment of the random measure N. The continuous 
quantum trajectory (p t ) satisfies 



Pt 



Po + 



+ 



L(p s .) +Tr[j(p s .)]p s ^ - J(p 

J(Ps-) 



ds 



J[0,K] 



Tr{J( Ps -)] 



Ps 



lo<x<Tr[j( PB -)]N(., ds, dx). 



(20) 



It is interesting to notice that we can work on [0, K] because Tr[J(p t )] < K for all process 
valued in the set of states. Based on this property, the following remark gives another 
equivalent way to define the solution of ©• 

Remark: The function t — > card(N(., [0, K] x [0,£])) = defines a standard Poisson 
process with intensity K. Thus for the filtration J-'t, we can choose the natural filtration 
of this process. The Poisson random measure and the previous process generate on [0,T] 
(for a fixed T) a sequence {(ri,£i),i G {1, . . . where each r, represents one jump 
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time of M. Moreover, the random variables £j are random uniform variables on [0, K\. 
Consequently, we can write our continuous quantum trajectory in the following way 



Pt 




The next section concerns the final convergence result. 

3 Approximation and convergence theorems 

This section is devoted to convergence theorems. The discrete quantum trajectory defined 
in Section 1 is shown to converge to the continuous quantum trajectory which is solution 
of the jump Belavkin equation ([20]) . 

The way to proceed to show this result is the following. Let h = 1/n be the time 
of interaction between Ho and a copy of Ti. This way, we shall define discrete quantum 
trajectories depending on parameter n. Next, based on result of Attal Pautrat in [3], 
we introduce asymptotic assumptions on the unitary operator U describing the evolution. 
Next, for special observable, we show that the corresponding discrete quantum trajectory 
is a "good candidate" to converge to the solution of the jump-equation. To prove the final 
result, it needs two step. 

• The first step is to realize the discrete quantum trajectory (in asymptotic form) in 
the same probability space of the continuous quantum trajectory. 

• The second step is to compare the discrete quantum trajectory with an Euler scheme 
of the jump-equation. 

Remark: In the literature, there exist classical Euler scheme results for jump-equation 
(see [7]). At our knowledge, the jump-Belavkin equation do not come into these classical 
setups. As a consequence, we establish the convergence of Euler scheme in our context. 
Let start by investigating the asymptotic behavior of discrete quantum trajectories. 

3.1 Asymptotic form of discrete jump-Belavkin process 

In this section, we present the asymptotic assumption concerning the interaction. Next we 
translate this assumption in equations 



C (p k ) + Ci(pk) + 



Qk+i 



£o(Pk) + 



Pk+i 



£i(pfc) X k 



(21) 



Pk+l 



Pk+l 



Qk+l 
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Consider a partition of [0, T] in subintervals of equal size 1/n. The unitary evolution 
depends naturally on n and we put 



U(n) 



L 00 (n) Loi(n) 
Lio(n) Lu(n) 



In the quantum repeated interactions setup [3], it is shown that the coefficients must 
obey precise asymptotic to obtain non-trivial limit when n goes to infinity. More pre- 
cisely, they have shown that quantum stochastic differential equations (also called "Hudson- 
Parthasarathy equations") , describing continuous time interaction models, can be obtained 
as continuous limit models of quantum repeated interactions by rescaling discrete interac- 
tions (such results use Toy Fock space and Fock space formalism). 

We do not need the total result of [3], we just need the expression of L^. In our context 
the asymptotic are the following 



I +-{-iH -)-CC* 
n 2 



o | - 

n 



Loo(n) 
L 10 (n) 

Concerning the description of the total Hamiltonian H tot (n), we can write 

I \ I 



L C + ofI 

n \n 



(22) 
(23) 



Htot(n) 



H ® I + I <g> 







c 





1 



I 




o I - 

n 



With the time discretization, we obtain a discrete process depending on n 
Pk+i{n) = £o(n)(pk(ri)) + Cx{n){p k {n)) 
+ 



C (n)(p k (n)) + 



Pk+i{n) 

Remind that by definition, for (X^), we have 



Qh+i(n) 



£i(n)(p k (n)) 



X k+1 (n). (24) 



X k +i(n)(i) 



Qk+i(n) 
Pk+i(n ) 

Pk+i(n) 



with probability p k+ i(n) if i = 
with probability q k+ x(n) if i = 1 



(25) 



where p k+ x(n) = Tr I®P U(n)(p k ® /3)U(n)*I (8)P = 1 —Qk+i( n ) (remind that P is one 
of the eigen-projector of the measured observable). Hence, the continuous limit behavior 
of the sequence (p k ) will depend on limit behavior of the random variables (X k ) and then 
depends on asymptotic behavior of probabilities p k +i and g^+i- 

By applying the asymptotic assumption in equation (l2"4l) . depending on the observable, 
we get two different situations. 
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If the observable A is diagonal in the orthonormal basis (Q, X) of H, that is A = 
1 0\ (0 

o o +Al V o i 



( n a ) + ^i ( n ? )• We obtain the asymptotic for the probabilities 



p k+1 (n) = l--Tr[J(p k (n))]+o(-) 
q k+ i(n) = ^Tr[J{p k {n))]+o(^j. 



The discrete equation becomes 



Pk+i(n) ~ Pk{n) = -L(p k (n)) + o(- 

71/ \ Th 



+ 



J{pk{n)) ( , ' 

-Pfc(n) + o(l) 



Tr(J(p k (n))) 



Vlk+i(n)p k+1 (n) X k+1 (n). 



Poo Poi 

PlO Pll 



• If the observable is non diagonal in the basis (Q,X), we consider P 

and P 1= ( %0 Qm ) we have 
V Qw qn J 

Pk+i = Poo + ^Trip k (p 01 C + p 1 oC*)] + -Trip k (po (C + C*))] + o(-) 

q k+ i = qoo + -^Tr[p k (q 01 C + q w C*)] + -Tr[p k {q m {C + C*))] + o ( - 

n \n 

The discrete equation becomes then 

Pk+i ~Pk = -L(p k ) +o(-) + [e l6 Cp k + e- i9 p k C* 

-Tr[p k (e* e C + e-"C*)] p k + o(l)] -±=X k+1 . 

From these descriptions, we can define processes p[ nt ] (n) by 

[nt]-l 

P[nt](n) = po+ ^ Pk+i - Pk 

k=l 

[nt]-l . [nt]-l 



Po+J2 ( L (^) + °( 1 ))-+ Qi(Pk)X k +i(n), 

k=l ^ k=l 



(26) 

where the expression of Q, depends on the expression of the observable (diagonal or not). 
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In the non diagonal case in [22J, an essential result is the proof of the convergence 

[nt] 

W n {t) = ^Y j X k {n) W t (27) 

V fe=l 

where T> denotes the convergence in distribution and (Wt) is a standard Brownian motion. 
Next by using a theorem of convergence for stochastic integrals due to Kurtz and Protter 
(cf [H],[T9]), it was shown that discrete quantum trajectories converges to solution of 
diffusive Belavkin equations. The convergence of (W n (t)), in the non diagonal case, is 
proved independently of the convergence of (p[ n t])- Besides this convergence allows to 
conclude to the convergence of (p[ n t]) with Theorem of Kurtz and Protter (such result 
needs namely the convergence of the driving process before to consider the convergence of 
quantum trajectories). 

In the diagonal case, we expect that the discrete quantum trajectory converges to the 
solution of jump-equation. In this similar result for driving processes would be 

N n (t) = J2x k (n) -^N t - Tr[J( Ps _)]ds. (28) 

k=l J0 

Apparently such result would need first the result of convergence of (p[ n t\) to the solution of 
jump-equation. Actually in the same spirit of definition of solution for the jump-equation 
involving the simultaneous existence of (N t ) and (pt) because of the dependance of the two 
processes, in order to prove the final convergence result, we have to prove the simultaneous 
convergence of (N n (t)) and (p\ n t])- Since we cannot show the convergence of N n (t), the 
convergence result of Kurtz and Protter cannot be applied. This justifies the use of random 
coupling method in order to compare in a simultaneous way all the processes. 

Before to deal with this method, we present the result of the convergence of Euler 
scheme. 



3.2 Euler-scheme for jump-Belavkin equation 



The literature abounds in references about Euler scheme approximation for stochastic 
differential equations (cf [7] , [12] , [T3] ) . The non-usual case of jump-Belavkin equations is 
not really treated, that is why we present the different result in this situation. 

When we want to study Euler scheme approximation for stochastic differential equa- 
tions, an important property is the Lipschitz character of the coefficients. Remind that 
our equation is of the following form 



p t = Po + f(p s ~)ds 
Jo 

J(Ps 



J[0,K) 



Tr[J(p s -) 



Ps 



l 0<x<Tr 



U(jt s _)]N(.,dx,ds) 



(29) 
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We must transform this equation to have Lipschitz property. We want to write it in 
the following way 



Pt = Po + f(n a -)ds+ / [q{n*-)] lo<x<TrU(n.-)]N(., dx, ds) (30) 

JO JO J[0,K] 

where / and q are Lipschitz functions and defined for all matrices. The consideration 
about the Lipschitz property is a pure technical aspect and can be admit by the reader. 
The Euler scheme is given by the formula ( l38l) below. 

Concerning /, we have seen that the solution of ( 1291) is obtained by truncation method 
because / is not Lipschitz but C°°. It was shown that the truncation is unnecessary since 
the solution is a process valued in the set of states. As a consequence we can consider that 
the function is truncated and then Lipschitz. We denote by F its Lipschitz coefficient. 

Concerning q, we must control the function defined on the states by 

J(p) 



9 ■ P 



P 



lo<Tr[J(j>)] (31) 



Tr[J{p)\ 

We transform the expression and define a function q which is C°° and such that 

g(p) = <?(p)lo<Tr[J(p)]- 

To construct the function q, it depends on the invertible character of C. 

If C is invertible, the function defined on the set of states by p — > Tr[J(p)} is continuous. 
With the fact that for all state p, we have Tr[J(p)} > and a compactness argument, the 
function p — ► J^(p)/Tr[J'(p)] is extendible by a function C°° defined for all matrices. 

If C is not invertible, there exists a unitary-operator V and two complex scalars a et 
j3 such that 

VCV*=(* (32) 

Before to go further, we have to show that the jump-equation is equivalent under unitary 
modification. For this, define for any unitary-operator V 

Jv(p) = VCV*(p){VCV*Y (33) 

f v (p) = -i[vHv*,p]-~{vcv*(vcv*y, P } 

+VCV*p(VCV*Y - Jv(p) + Tr[J v (p)]p (34) 
Jv{p) , ori 

*-0<Tr[J V (p)]- 1<50J 



9v{p) 



[Tr[J v (p)} 

We have the following proposition which expresses unitary equivalence. 

Proposition 4 Let V be any unitary operator and let (p t ) be the solution of the jump 
Belavkin equation, then the process {^t '■= Vp-tV*) valued on the set of states satisfies 

7t = 7o+ / fv{ls-)ds+ / / [gvils-)) l <x<Tr[j v ( ls -)]N{.,dx,ds). (36) 

</0 ^0 J[0,K] 
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The proof is a straightforward computation. Such unitary equivalence allows us to 
transform g without change Lipschitz property of /. 

Now let see how we can construct the function q in the case where C is not invertible. 
Let V be the unitary operator involved in expression (1321) . we get 



9v(p) 



1 




P 



l 0<Tr[J v (p)}- 



Hence, the expression of q is clear and by using the unitary transformation given by Propo- 
sition 4, we can consider the following equation 

Pt = Po+ f(ji a -)ds+ / [qins-^lo^TrU^^Nt^dx^s) (37) 
Jo Jo J [0,1] 

which admits a unique solution by Theorem 4. 

Another time, the fact that q is C°° allows to consider that it is a Lipschitz function 
with a truncature method. We denote by Q the Lipschitz coefficient of q. 

With these technical precautions, we can consider the Euler scheme 

fc+i 

9 k +i = k + -f(6 k ) + / [q(9 k )]l < x <Re(Tr[j(e k )])N(.,dx,ds). (38) 

n J± J[0,K] 

Let fixe an interval [0, T] and for all t < T, we define k t = maxjfc 6 {0, 1, . . .}/- < t}. For 
all t in ]|, £±i], we put 

t = k + [ f(9 k )ds + [ [ [q{9 k )]l o < x < Re{Tr[J{0k)]) N{.,dx,ds). (39) 
<M 7 a 7 [o,i] 

It is worth noticing that we have 9k = 9 k for all k. We then have for t < T 

n 

fct-l ,.i=±I fet-1 ,.*±1 

'[0,K] 



•VI J- /• 1 rvt J- /» 1 — /» 

6t{n) = Y] / f(9 k )ds+y2 / [g(6 ) fc)]lo<x<i? e (Tr[ l 7(e fc )])iV(-,^,rfs) 

fc=0^ fc=0^ ■'P.-Kl 

+ / f(9 kt )ds+ / / [g(4J]lo<x<i?,e(Tr[j(e fct )])A^(.,rfx,rfs). 

At •-'fct ^[0,K] 

Likewise for the solution (/i t ) of the Belavkin equation, we can write 

fct-l »i=±I fct-l 



V / f(p s -)ds + V / / [g(/i s _)]l < a; <Tr[j(^_)]^(.,rfx,rfs) 

/ f{n,-)ds+ / / [g(/i s _)]l <x<Tr[J( Ms _)]A^(.,rfx,c/s). (40) 
At Jk t J\0,K] 



' [0,K] 

Before to express the convergence theorem, we need the following proposition. 



21 



Proposition 5 Let (fi t ) be the solution of the jump-Belavkin equation, then there exists a 
constant M such that for all (s, t) G 



E[||/i t -//J] < M\t-s\. 



(41) 



Proof: From the fact that the solution of the jump-Belavkin equation is valued on the set 
of states, we have for all t > 0, the quantity \\fi t \\ < 1 almost surely (we do not specify the 
norm because we just need the fact that the solution is bounded). For < s < t, we have 



IH-Va= \ f{li u -)du+ / / [q{ii u J)\lQ< x < Tr [j^_)]N(. : dx : du). 

Js Js J [0,1] 



(42) 



By using the property of random Poisson measure, in particular the property about the 
intensity measure we have for < s < t 



BIW^-^W) < E 
+E 

< 



f({i u J)du 



s J[0,K] 

E[||/(^_)||]d« 



b(Vu-)}lo<x<Tr[J(^)]N(., dx, du) 



+E 



< 



< 



s J[0,K 

E[\\f(^)\\]du + E 



[q(^u-)}lo<x<Tr[j(^ u .)]\\ N(., dx, du) 

/ / \\q(Hu-)lo<x<Tr[J(tM 
Js J\0.K] 



dxdu 



sup \\f(R)\\+2)du 

\\m<i 

< M(t-s) 
where M is a constant. The result is then proved. 



□ 



We can now express the theorem concerning the convergence of the Euler scheme. The 
particularity of the jump-equation comes from the stochastic intensity depending on the 
solution. 

Theorem 3 Let T > 0, let (9t) be the process (Iff 9f\ constructed by the Euler-scheme on 
[0,T], let (fit) be the unique solution of the jump-Belavkin equation §Efi§ . 
We define for u < T and n large enough 



n 



E 



sup 

0<t<u 



6 t (n) - fi t 



(43) 



22 



So there exists a constant T which is independent of n such that for all u < T 



ZJn) < 



n 



(44) 



Let X>([0, T]) denotes the space of cddlag matrices processes endowed with the Skorohod 
topology. 

Hence the Euler scheme approximation (9 t ) converges in distribution in X>([0, T]) for 
all T to the process-solution [fi t ) of the jump-Belavkin equation. 

Before to give the proof, it is interesting to compare this result with the classic ones. 
In the literature, similar result appears often by using a L2 norm [7], that is terms like 



E 






2" 


sup 


6 t (n) - fi t 






_0<t<u 







are usually considered (to apply Ito Isometry result). Results like Z u {n) < Y/n 2 are 
then shown and almost surely convergence are then obtained. In the following proof, it is 
explained why such norm is not appropriate in our case. 

Proof: The equations concerning the Euler scheme and the solution of the jump-Belavkin 
equation give us the following formula 



9 tin) - fi t 



k t -i fc+1 t 
E / n [f(0 k )~ f(fi s -)]ds + I [f(9 kt )-f(^)]ds 

k = n 

kt-1 »fc±l 

+ ZJ / / [q(8k)lo<x<Re(Tr[J(9 k )]) 



kt 



k=0 



[0,K] 



+ 



kt J{0,K] 



-q(ti>s-)lo<x<TrlJ{n a -)])N(., dx, ds) 

(l{6kt)1-0<x<Re{Tr[J(e kt )\) 

-g(/i s _)l < a; <Tr[j( /t3 _)])^(., dx, ds). 



Let consider for u < T, the quantity Z u {n) 



E 



u 



9 t {n) - fi t 



Let consider 



separately the drift term and the term concerning the random measure. For the drift term, 
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by the fact that / is Lipschitz, we have 

k t -i n h+l 



E 



< E 



sup 

0<t<u ^ 
'k u -l „fc±l 



\\f(O k )-f^ s -)\\ds+ / \\f(O kt )-f(pi s _)\\ds 



E 



. k=0 n 
k=0 J n 

k u -l n h+l 



kt 



\f(O k )-f(^.)\\ds + / /(^)-/(/x s _) 



fc=0 

fc„-i -.Mi 



* E 

k=0 

k u 

+ E 

fc: 

ku~ 1 

* E 

fc=0 

+ E 

h 

ku 1 

s E 

k=0 

ku 

+ E 



E 



f(9 k )-f(^)\\]ds+ f e[ fiO^)- f{^)\\ds 

J K u 

fit-is-) - /(/•***) ds 



FE 



fc=0 

k u -l „h+i 



k=0 



E 



f(»s-) - 


f(Vk)\ 

n | 


\ ds+ L 


k_ — jd.k 


ds + 


/ FE 


n n 




ku 




n 


ds + 
J Jk 



E 



ds 











pu 








L " FE 


sup 


O t - Ht 




ds+ FE 


sup 


Ot ~ Ht 




n 


_0<t<s 






J k u 


_0<t<s 







k=0 

ku-l -M±l 



k u — /ifc 

1 

f(Vs-) ~ ) 

n 

ds 



E 



ds 



k=0 



FM [ s --)ds+ I FM is- — )ds 



n 



J ku 



n 



< A 



1 



Z s ds H — ) (A is a suitable constant). 



The analysis of the random measure terms is more complicated. Let fixe an index k, 
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thanks to the properties of random measure, we have 



E 



< E 



< E 



fc+1 
n 

k '[0,K] 

fc+i 



[0,K] 



+E 



fc+1 
n 

I [0,K] 



< E 



k >[0,K 

fc+i 



+2E 



[0,K] 

fc+i 



[q(0k)]lo<x<Re{Tr{J(8 k )]) ~ [<?(/^- )] ^0<x<Tr[ J{n„- )] N(. , d,X, ds) 

[q(6k)}~lo<x<Re(Tr[j(e k )\) — [q(fJ, s -)] lo<x<Tr[j(^ s - )] 1 1 N(.,dx,ds) 

[q{&k)\^Q<x<Rs{Tr[J{e k )\) ~ [<l(V>s-)] ^-0<x<Tr[J(9 k )} || N(.,dx,ds) 
l|[<?(^s-)]il X \la<x<TrU(n a -)] - lo<x<Re(Tr[J(8 k )})\ N(.,dx,ds) 

q(9 k) — q(fi s _) N(.,dx,ds) 
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\^-Q<x<Tr[J(jjt a -)] — lo<x<Tr[J(6 k )}\ N(.,dx,ds) 



>[0,K] 

As q is bounded by 2 on the set of states, we have 



(45) 





fc+1 




M < E 


r/ « 











-2E 



fc+i 



[0,X] 



iV(., efcc, o?s) 



L O<a:<max(Tr[J(/i s _)],i?e(Tr[ 1 7(0 iL )])) 



1 Q< !r <inin(Tr[J'( i u s _)],-Re(Tr[ l 7(0 A )])) ) N Vi dx ' ds ) 



Furthermore, we have Re(Tr[J'(ji B -)]) = Tr[J~(ji s -)] for all s. Hence, by linearity and 
continuity of the function trace, for any matrices A and B, there exists a constant R such 
that 

\Re(Tr[J(A)]) - Re(Tr[J(B)})\ <R\\A- B\\ . 
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It implies 





" r h+l 








r r^±k 


< E 


h Q 

^ n 


9 k — A* s - 

n 


ds 


+ E 


I: 



Re(Tr[J(9 k )])-Tr[J{n s -)) 



ds 



< 



< 



(R + Q)E h -fx 



fc+i 



(R + Q)V 

< (R + Q) 



fc+i 



Z s ds + — 
n z 



ds 
ds 



fc+i 



— Ll s - 



(R + Q)E 
we do the same as the drift term). 



(46) 



The term between k u and u can be treated in the same way. By summing, we obtain 
finally the same type of inequality for the term with the random measure. As a consequence 
there exist two constants Fi and F 2 which depend only on T such that 



Z U <F X Z s ds + 



n 



(47) 



The Gronwall Lemma implies that there exists a constant V such that for all u < T 

r 

Z u (n) < - 
n 



(4f 



where T is a constant independent of n. The convergence in 2?([0, T]) is an easy conse- 
quence of the above inequality. The result is then proved. 

Now we can justify why the L 2 norm is not appropriate to deal with such equation. 
Indeed, the last term of (145p involves an integral with a difference of two indicator functions. 
The difference of the indicator functions is equal to zero or one and it gives the same result 
if we choose the L 2 norm, that is 



^Q<3s<Tr[J{fs-)] 



\-0<x<Tr[J(e k )] | 



<-0<x<Tr[J(fj, s -)] 



tO<x<Tr[J(e k )\ | 



After in the above proof, the integral of this term is calculated and if we use L 2 norm, we 
loose the homogeneity in term of L 2 norm. As the final result relies on Gronwall lemma 
lemma, this homogeneity is actually necessary to obtain an appropriate estimation. □ 

In the following section, we compare the discrete process with the Euler scheme. 
3.3 Convergence of the discrete process 

This section is devoted to the random coupling method of the discrete quantum trajectory 
and the continuous quantum trajectory. 
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Consider the probability space (fl, T , F t , P) where the solution (p t ) of the jump Belavkin 
equation (I2U1) is defined. Let construct the discrete quantum trajectory in this space. Let 
n be fixed, we define the following sequence of random variable which are defined on the 
set of states 

fk+lfaw) = ljV(w,Gi(f,))>0 (49) 

where G k ( V ) = {(t, u)/ £ < t < *±I, < u < -nln(Tr[jC (n)(v)})}. 

Let p — P be any state, we define the process (p*,) for k < [nT] by the recursive formula 

C (p k ) + C x (p k ) 



+ 



+ 



(%fi(pfc, •) - Tr[£i(p fc )]) . 



(50) 



Tr[C (p k )} Tr[C 1 (p k )}_ 

This random sequence and the operators Ci(n) depend naturally on n following the asymp- 
totic of the unitary evolution. Thanks to the Poisson distribution property, the following 
proposition is obvious. 

Proposition 6 Let T > be fixed. The discrete process {p k ) k <,[ n T} defined by ( \50§ has 
the same distribution of the discrete quantum trajectory (p k ) k <[ n T} defined by the quantum 
repeated measurements principle. 

This proposition is a consequence of the fact that for all Borel subset B G £>(R 2 ), we 

get 

A(B) k 



P[N(B) = k] 



k\ 



exp(-A(5)), 



where A denotes the Lebesgue measure. 

In (fl, J 7 , jF t , P), the process (p k ) satisfies the same asymptotic than the discrete quan- 
tum trajectory, that is, 



Pk+l — Pk 



Pk 



j{pk) 



Pfc + °p fc (l) 



Vk+l(Pk>-)- ( 51 ) 



[Tr(J(p k )) 

Before to compare the discrete process fl50|) . we need another process which concerns the 
approximation of the intensity. In (fl, T, T%, P), we define the random variable sequence 
defined on the set of states 

v k+1 (r],u) = lN(b>,H k ( v ))>o (52) 

where H h (rj) = {(t,u)/% <t< *±1,0 < u < Tr[J(rf)]}. Let p = p be a state, we define 
the following process in (fl, T , J-'t, P), for k < [nT], we put 

£o(p~k) + £i(p fe ) 

£o(Pfc) A(p fc ) 



Pk+l 



+ 



+ 



Tr[C (p k )] Tr[A(p*)] 
Hence, we have the same asymptotic form 



("k+i(p 



ki 



Tr[dM]). 



(53) 



Pk+i ~ Pk = -if(Pk) + °P fe (l)] + 



n 



J{pk) 
Tr(J(p k )) 



Pk + °P k iX) 



"k+i(p. 



fej 



(54) 



Regarding the process (p k )o< k <[ n T] and (p k )o<k<[nT\, we have the following proposition. 
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Proposition 7 Let (pk)o<k<[nT] be the discrete quantum trajectory defined by the formula 
( \51\ and let [Pk)o<k<[nT] be the sequence defined by the formula (\5'J[\ . Let assume that the 
two sequences are defined by the same initial state p. 
For k < [nT], we define 



AAn) = E 



We have for all k < [nT] 



sup \\pi{n) -Pi{n) 

0<i<k 



A k (n) < o - 
n 



where the little o is uniform in k. 

Proof: Remind that the discrete quantum trajectory satisfies 



pk+l ~~Pk = ~[f(Pk) + °p h 0)] + 



J(Pk) 



h+i(p~k,-)- (55) 



[Tr[J(p k )\ 

We can remark that all the rest °pfc(l) are uniform in k because the process (pk) is valued 
on the set of states and so is bounded. Hence, we can write this equation in the following 
way using / and q 



Pk+l ~ Pk = ~[f(Pk) + °p fc (l)] + [q{Pk) + Op fc (l)]z> fc+1 (p fc , .)• 



(56) 



We have a similar asymptotic form for the process f[5ol) . As a consequence, we can compare 
the two processes 



Pi- Pi 



n 



Hence, we have 



SUp ||p; - Pi\ 

0<i<k 



i-1 

£ 

3=0 
i-1 



fc-i 1 

^ E^||(M-)-/(Pi)+o^(i) 

j=0 
k-1 

+ Yl WWi) + °^(i))^+i(Pi, •) - q{PjPj+i{Pj, •) 
fc-i 



< 



j=0 

fc-i - 

j=0 j=0 

fe-i 



(57) 
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By defining the filtration Q 3 = cr{i>k(pk-i, .), Vk{Pk-ii •)■> < ^ — j} f° r 3 > 0; we have by 
the independence of the increments of a Poisson process 



E 

= E 
= E 

= E 



\\q(f?) + 0^(1)) - q(p 3 ) - o- pj (l)))V j+1 (p j: 



E [Hp,) + o^.(l)) - q {- Pj ) - o- Pj {i)))v j+1 {- Pj , .)\\/g 3 

]\q(p 3 ) + o^.(l)) - q{- pj ) - op (1)))E Oil/ft-] 



||g(p,) + o^.(l) - q{- pj ) - 0^ (1)|| 1 - exp --Tr^fo)] 



< E[Q||p,-pJ](irr[^01+oQ))+o(i 



because all the rest are uniform in j. The same way, by using the filtration, we get 

E \\(Q(Pj) + •) - -)| 

= E[ 11^) + o^.(l)) || E [\i> 3+1 (p 3 , .) - .)| /$] 

For the second term of the product, by definition, we have 

E \vj+i(Pj,-) -Vj+i(Pj,-)\ /Gj 
= E l N (., Gj fa))>o ~ 1jv(.,^(^))>o| 

= E l{N(.,G j (p j ))>0}A{N(.,H j (p j ))>0}/Gj\ 

= P [{N(., G 3 (p 3 )) > 0} A {N(., H 3 {- Pj )) > 0}/Q 3 ] . 



(58) 



We denote by W 3 = {(t,u)/i < t < i±I, min (Tr[J(p 3 )\, -n\n(Tr[C (p 3 )])) < u < 
max (Tr\J{j) 3 )\, — n\n(Tr[£ (pj)]))}. Hence we have 



E 



= P[i Wj> olGi\ 

= 1 - exp ^ - i(min {Tr[J(p 3 )}, -n\n(Tr[£ (p 3 )})) 

,HTr[£ (p 3 )})) ) 



max . . . , — n J 



1 — exp 

n 



Tr[J(pt)}+nln[Tr[£ (p 3 )} 



n 



Tr[J(p 3 )]+n\n(Tr[C (p 3 )}) 



I - 

n 
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Besides we have Tr[£ (p~j)] = pj + i = 1 — -Tr\J{pj)] + o(-), hence 

1 



E 



n 



TriJip^-TriJfa)) 



+ o 



(59) 



As (pk) is a process valued in the set of state, it is uniformly bounded, we then have 



E 



\(q(pj) + o Pi(i))(^i+i(Pi, •) - 













< K~E 


\\pj-pj\\ 


-a 



(60) 



By taking expectation in ( 1571) and using the inequalities ( |58ll 591) . we obtain finally the 
inequality 



fc-i 



\Pj-Pj\ 



^— ' n i 

j=0 



+ o 



(61) 

The result follows with a discrete Gronwall Lemma. □ 
Now, we can compare the process obtained by the Euler scheme and the process defined 
by the formula floT|) . The result is resumed in the following proposition. 

Proposition 8 Let (j> k )o<k<[ n T\ be the process defined by the formula ( [5^] ) and let (0k)o<k<{nT] 
be the process obtained by the Euler scheme of the jump-Belavkin equation. Let assume 
that the two sequences are defined by the same initial state p. 
For k < [nT], we define 



S k (n) 



E 



sup \\9i(n) - p^n) 

0<i<k 



Hence there exists a constant F such that for all k < [nT] 

F 

S k (n) < —■ 
n 

The proof is based on the Gronwall Lemma but it uses finer property of the random 
measure induced by the Poisson point process. This is a generalization of the Poisson 
approximation studied by Brown in [6]. 

Proof: Thanks to the fact that random sequence (pj) is bounded, the o-p. = °(1). It 
implies that for i < k < [nT] 



0i ~ Pi 



E^[M-)+°(i)-/(Pi)]+E ' 

i-1 
3=0 



[0,1] 



[g(6' j )]l < x <Tr[j(e J )iV(., dx, ds) 
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We treat the random measure part and the drift term part separately. Let us denote 
S k = E [supo^fc pi - p^W , we have 



S k < E 



E^I|[M')-M-) + o(l)]|| 

j+l 

I I Q.^j)^o<x<Tr[j{fi j )\N{., dx, ds) - (q(pj) + o(l)) V j+1 (p j: .) 



k-l 

E 

.3=0 



As / is Lipschitz, we have 



E 



E^llM)-M)+°(i)ll <^E^+°QV 



For the second term, we get 



E 



E 



'ifc-i 

E 

U=o 
fe-i 



2+1 



[0,1] 



< E 



+E 



fc-i 



E Ito w> ^to)) - to) + ° (!)) • 

fe-i 

E lk(^)*W^> •) - + ° (*)) 77 i+i(Pj> 

j=0 
k-l 



.3=0 



< ^E[(g||^-p-|| + o(i))x|i7 j+1 (p.,.)|] 

3=0 
k-l 

+ E E [II^)II x IN^HM))-*^,.)]]. 

3=0 

Here we introduce the following discrete filtration 

^ = ffWft-i,-),iv(..W<i} 



(62) 



It allows to compute the previous terms. It is clear that the random variables p • and 9j 
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are Tj measurable, we then have 



K — l 

E - fate) + (!)) 

.j=o 
fc-i 

< E E [Wll^ " + (!)) >< E [fc-nfo, 
fc-i 

+ E E 01^)11 x E[|7V(.,tf j( 0.)) -I7, +1 (p.,.)/^|]] 



By conditioning 
Hence we have 



with respect to J-'j, the random variable Vj +1 {p^.) is of Bernoulli type. 



E 



n J 



. n 



For the second part, we have almost surely 

E[|iV(.,^(^))-^ + i(P,,-)/^l] 

< E [{N^H^)) -N^Hfc))/?,] +E [N^Hjfa)) -V j+1 (p j: .)/^] 

< \ \Tr[J(p 3 ))} -Tr[J(ej)]\ +E [N^Hjfc)) -V J+1 (p 3 , .)/^} 
A u n „ A' f l\ 

< -\\pj-0j + — + ° — ■ 

The o are uniform in j because (Pj)j is uniformly bounded. For the second term, the 
above inequalities and the fact that the Euler scheme is bounded implies that there exist 
two constants K\ and K 2 such that 

k-l 

[11^)11 x E^H^)) -^ + i(Pp-)\ m] 

3=0 

< K £l Sj + K> + a (L). 

3=0 v 7 

For the first part, we have an equivalent inequality. Thus we can conclude that there exist 
two constants G\ and G 2 such that 



s k < Gl ! £is j + ^ + o(l). 

3=0 v 7 



(63) 
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The discrete Gronwall Lemma implies that there exists a constant F independent of n such 
that for all k < [nT] 

S k < -. 

n 

The proposition is then proved. □ 

By using this two properties we can now express the final theorem. 

Theorem 4 Let T > be a fixed time and let (Q, J 7 , Ft, P) be the probability space of the 
poisson point process N . Let n be an integer and let (p[nt])o<t<T be the discrete quantum 
trajectory defined for k < [nT] by the equation 



Pk+l 



C {p k ) +Ci(p k ) 
£o(Pfc) 



+ 



+ 



Tr[C (Pk)] TrlC^W 



{vk+i{p~k, •) - Tr[Ci(p k )]) . 



Let (pt)o<t<T be the quantum trajectory solution of the jump Belavkin equation on [0, T] 
which satisfies 



I't 



t 

//,:, ' / ./'(/',. )r/.s 

'o 

" ' r J(Ps-) 



+ 



J[Q,K] 



[Tr[J(p s -)} 



Ps 



1-0<x<Tr[J(Vs-)] N (-> ds). 



If Po = Po> then the discrete quantum trajectory (p[ n t\)o<t<T converges in distribution to 
the continuous quantum trajectory (pt)o<t<T i> n D ([O,? 1 ]) for all T. 

Proof: Let n be large enough. For k < [nT], we define Pk = Ph.- We define 



B, = E 



sup 1 1 pi 

0<i<fc 



Pi 



Thanks to Proposition 8 and Theorem 5 concerning the Euler scheme, there exists a con- 
stant R independent of n such that for all k < [nT] 



R 

B k < -. 

n 



(64) 



It is worth noticing that the process (p[nt])o<t<T converges in distribution to (pt)o<t<T 
for all T in X>([0, T]). Thanks to this fact and the inequality the convergence in 

distribution of (p[ n t})o<t<T to (p t )o< t <T is proved. □ 
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